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Abstract: We revisit the problem of why stars become red giants. We modify the physics of a 
standard stellar evolution code in order to determine what does and what does not contribute to a 
star becoming a red giant. In particular, we have run tests to try to separate the effects of changes in 
the mean molecular weight and in the energy generation. The implications for why stars become red 
giants are discussed. We find that while a change in the mean molecular weight is necessary (but not 
sufficient) for a IMq star to become a red giant, this is not the case in a star of 5M©. It therefore 
seems that there may be more than one way to make a giant. 
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1 Introduction 

The question 'why do stars become red giants?' is per- 
haps one of the longest standi ng problems in stellar 
astrop hysics. In a recent paper ISugimoto fc Fuiimotol 
(2000) provided a long list of work on the subject, with 
publication dates spanning over four decades. The 
problem has been approached from many different an- 
gles, from considerations of polytropic solutions (e.g. 
lEggleton fc Cannonlll99ll ) to detailed numerical mod- 
elling (e.g llbenll993f l. Despite all the investigation into 
the subject, the question has yet to receive an answer 
that is satisfyingly simple and sufficiently rigouroufl 
There is still no consensus on why stars become red 
giants. Theories to explain the phenomenon are many 
and varied. Some assert that a 's oftening' of the effec- 
tive equation of state is re quired {Egglcton fc Cannon! 
ll99ll ; lEggleton" et al .1 19981 ) . others suggest that a strength- 
ening of the centra l gravitational field is required (jHoppner 
11973 : IWeis j| 1983T ) . Some attribute the behaviour to a 
gravothermal instability in the core ljlbenll993h . others 
suggest it is a therma l instability in the stellar envelope 
(R enzini et al.lll992T ). The list of potential explana- 
tions is long and to review them all would take more 
pages than t he authors have been permitt ed for this 
contribution. ISugimoto fc Fuiimo to (2000) give a re- 
view (as well as a rather fierce rebuttal) of the leading 
theories and the interested reader is directed to this 
paper for a more thorough discussion. 

Here we shall only briefly review some of the work 
that has been don e on the problem via direc t numer- 
ical experiments. iHoppner fc Wei gcrt (1973]) investi- 
gated the effects of an increased gravitational field 
in models of 1.4 and 1.8 Mq, as well as looking at 
the effects of inhomogeneities in the mean molecu- 
lar weight and in energy generation. They concluded 



x Of course, we are labouring under the assumption that 
such a solution does indeed exist. 



that only a strong gravitational field could p roduce 
prope rties similar to red giants. Subsequently, IWeissI 
(|1983| ) extended their work to cover stars in the mass 
range 1 < M/Mp ) < 8 . Using polytropic models, 
iFrost fc Lattanziol l| 19921 ) later demonstrated that this 
could not be the sole cause. 

IRenzini et all (|l992h suggested that stars become 
red giants because of a thermal instability in their en- 
velopes. In their view, expansion is initially driven by 
the envelope maintaining thermal equilibrium in re- 
sponse to increasing luminosity from the core. This 
expansion leads to local cooling and the recombina- 
tion of heavy elements. An increase in the opacity 
results, trapping energy in the envelope. This leads 
to a runaway expansion th at br i ngs th e star to the 
red giant branch. However, llbenl (|1993| ) subsequently 
showed that the opacity could not be responsible for 
the trans ition to a red giant structure. He did this by 
4^Hg models in which the opacity was held con- 
stant throughout the star. These models still became 
giants. 

In reviewing the literature, what is perhaps most 
striking is how contradictory much of the work is. For 
example, numerical experiments bv llbenl ljl993h show 
tha t the opacity does not p l ay a k ey role in the process, 
yet ISugimoto fc Fuiimotol (|2000h assert that "A key 
role is played by the gradients of opacity" . One of the 
reasons for this may be that most studies are limited 
in their scope, looking only at models of a particular 
mass. Studies examining a range of masses are the 
exception, rather than the rule. A systematic study is 
clearly warranted. 

In this work, we adopt the commonly-used ap- 
proach of running experiments using a detailed evolu- 
tion code. By altering the input physics in a controlled 
way and seeing what effect this has on the evolution, 
we hope to gain some understanding of what factors 
control whether a star becomes a red giant or not. We 
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emphasise that by the phrase "becomes a red giant" 
we mean that the star develops a condensed core with 
an extended envelope, moving to cooler effective tem- 
peratures as it does so. 

2 The stellar evolution code 

Our calculations have been carried out using the STARS 
stellar evolution code originally developed by Egglcton 
(1971) and updated by many authors (e.g. Pols et al. 
[1993). The code follows the evolution of seven ener- 
getically important species, namely 1 H, 3 He, 4 He, 12 C, 
14 N, ie O and 20 Ne. A detailed d escription of the code 
and its features can be found in IStancli ffc (2006) and 
references therein. 

All our experiments have been performed on a 1 Mq 
model of metallicity Z=0.02. In our standard run we 
do not consider the pp-chains. The reason for this 
is that in some of the experiments the use of the pp- 
chains leads to models that do not have giant-like com- 
position profiles: they are not shell-like. The hydrogen 
abundance declines slowly towards the interior, rather 
than having a sharp drop at a particular location. By 
using only the CNO-cycle reactions to burn hydrogen, 
we obtain much sharper, shell-like profiles. The use of 
only the CNO-cycle reactions presents us with a prob- 
lem: the more concentrated energy release leads to the 
formation of a convective core. As we want the star to 
behave in a similar way to a normal 1 Mq model, we 
prevent any mixing from occurring in the convective 
corfl Figure [1] shows a Hertzsprung- Russell (HR) di- 
agram for our standard, CNO-burning only model and 
a normal 1 Mq model that includes both the pp-chains 
and the CNO-cycle. 
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Figure 1: HR diagram showing the evolution of a 
regular 1 M© model (solid line) and our standard 
1 M Q model which does not include the pp-burning 
reactions. The standard model is hotter because 
the star must contract more before the CNO ele- 
ments settle to their equilibrium value. 



2.1 Suppression of surface convec- 
tion 

In this model, we suppress the occurrence of surface 
convection by forcing the code to use the radiative 
temperature gradient and preventing any mixing of the 
chemical elements. The resulting evolutionary track is 
displayed in Figure [5] The starting model was taken at 
the ZAMS of the standard model. The model makes a 
sudden transition to lower temperatures at the begin- 
ning of the sequence because the star has to recover 
from the sudden loss of its surface convection zone. 

Rather than evolving to higher luminosities after 
the end of the main sequence, this model proceeds to 
lower surface temperatures. The radiative envelope 
does not transport energy efficiently like a convective 
envelope does. Instead, the radiation is absorbed and 
the star's radius increases, pushing the star to lower 
surface temperatures. This star still becomes a giant, 
though it does not reach as high a luminosity as the 
standard model. The existence of surface convection 
is not important for a star becoming a red giant. 

2 It would also be desirable to force the code to use the ra- 
diative temperature gradient, essentially removing convec- 
tion completely from the core. Attempts to do this proved 
unsuccessful as the models failed to converge. 
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Figure 2: HR diagram showing the evolution of the 
model without surface convection (dashed line). 
The standard model is also displayed for reference 
(solid line). 
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Figure 3: HR diagram showing the evolution of the 
homogenous model (dashed line). The standard 
model is also displayed for reference (solid line). 

2.2 Homogenous evolution 

In this model, we force the whole star to be fully mixed 
throughout its evolution. The resulting evolutionary 
track is shown in Figure [3] In this model, the star 
does not become a giant - its radius remains around 
one solar radius until hydrogen has been exhausted 
and then it begins to fall. The star moves blueward, 
not redward, throughout its evolution. 

2.3 The role of the mean molecular 
weight 

The homogenous model suggests that the mean molec- 
ular weight may play a role in a star's journey to gi- 
anthood. To investigate this, we have made a set of 
models that do not convert hydrogen to helium in the 
normal way. Instead, we burn hydrogen into a false el- 
ement which we shall call pseudohelium. This pseudo- 
helium is treated like normal helium in all respects but 
one. While it contributes to the opacity in the same 
way as helium and has the same ionization states and 
number of electrons as helium, we assign it a differ- 
ent atomic mass. By running various model sequences 
with different atomic masses for the pseudohelium, we 
can examine the effect that the mean molecular weight 
has on whether a star becomes a red giant. 

We have run test cases where the mass of the pseu- 
dohelium is I, 1.5 and 2. The results are shown in Fig- 
urefl] It should be noted that the case M psou doHc = 1-5 
keeps the mean molecular weight of the star constant 
(assuming full ionization) because the pseudohelium 
has two electrons associated with it. If the atomic 
mass of the pseudohelium is below 1.5, the star does 
not become a giant nor does it become red. For a pseu- 
dohelium mass of 2, the star does evolve into a red 
giant with a compact core and an extended envelope. 
Experimenting with pseudohelium masses between 1.5 
and 2 suggests that this transition is smooth. 

This sequence of tests, together with the earlier 
homogenous model, suggests that the mean molecular 
weight has a key role to play. If there is an insufficient 



Figure 4: HR diagram showing the evolution of 
the models with pseudohelium of varying atomic 
masses. The standard case is displayed for com- 
parison and corresponds to M psGU doHc = 4. 

change in the mean molecular weight throughout the 
star, the star does not become a giant. We shall return 
to this point later. 

2.4 The role of energy generation 

We should not expect a change in the mean molecu- 
lar weight to be solely responsible for a star's transi- 
tion to giant status. This is evident when one consid- 
ers the transition a star makes from the giant branch 
to the horizontal branch. The burning of helium to 
carbon and oxygen raises the mean molecular weight 
(even more so than the conversion of hydrogen to he- 
lium) yet the envelope shrinks. The nuclear burning 
also increases the energy generation in the core. One 
may therefore ask the question "What role does energy 
generation play in a star's transition to a giant?" . To 
investigate this, we make the following tests. 

2.4.1 Composition-independent energy gen- 
eration 

In this sequence, we allow hydrogen to be converted 
to helium in the usual way. However, we do not allow 
this process to generate energy. Instead the energy 
generation is given by 

e = e(p,T,X H = 0.7) (1) 

i.e. we generate energy using the current temperature 
and density of the model, but assume that the hydro- 
gen abundance is the ZAMS value. The energy gen- 
eration is thus decoupled from the chemical evolution: 
hereafter we refer to this model as the 'energy decou- 
pled' case. The energy generation of the model remains 
centrally concetrated throughout the evolution. 

The evolution of this model is shown in Figure [5] 
This initially looks promising: the model moves red- 
ward and we find that its radius is indeed increasing. 
However, closer inspection of the model shows that we 
have not made a giant in the sense that we wish. If we 
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Figure 5: HR diagram showing the evolution of the 
energy decoupled model (dashed line). The stan- 
dard model is shown for comparison (solid line). 
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Figure 6: The energy generation (solid line) and 
neutrino loss (dashed line) profiles as a function of 
mass adopted for the 'Forced Profile' run. 



plot the radius profiles of this model, we find that all 
layers are expanding as the model evolves. This does 
not give a structure that has a compact core with an 
extended envelope - our definition of a giant. This star 
is more akin to a star on its pre-main sequence. 

2.4.2 Forced profile 

The energy decoupled model has one unintended fea- 
ture: the energy generation rate (and its profile as a 
function of mass) changes with time. It is desirable 
to separate out these changes from the chemical ones. 
To do this, we take the energy generation profile (and 
the neutrino loss rate) as a function of mass at the be- 
ginning of the main sequence and impose this profile 
throughout the rest of the evolution. These profiles 
are displayed in Figure [H] The chemical evolution is 
allowed to occur as the usual function of temperature, 
density and composition. 

The evolution of this model in the HR diagram 
is shown in the top panel of Figure [7] The model 
does not become a giant. It evolves to the blue, not 
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Figure 7: Top panel: HR diagram showing the 
evolutionary track of the model with the energy 
generation profile fixed to its ZAMS value (solid 
line). The standard model is shown for compari- 
son (dashed line). Crosses mark the points in the 
evolution at which the profiles in the lower panel 
are taken. Lower panel: Radius, temperature and 
density profiles for the model at the points denoted 
by crosses in the upper panel. 
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the red, with all layers of the star contracting. The 
model has been evolved up to the point where the inner 
0.7 Mq has become hydrogen-exhausted. The temper- 
ature, density and radius profiles for this model show a 
distinct core has been formed (Figure [7] lower panel). 
The increase in the surface luminosity of this model is 
due to the release of gravitational energy as the star 
contracts. 



3 Discussion 

Of the above tests, two things stand out. First, the ho- 
mogenous model and the pseudohelium model which 
keeps the mean molecular weight constant do not be- 
come giants. Both these models evolve blueward with 
little change in their radii. This is despite their hydrogen- 
burning luminosities increasing by over two orders of 
magnitude from the ZAMS value. We may conclude 
that the mean molecular weight plays some role in 
whether a star becomes a red giant or not. Secondly, 
the forced profile model also does not become a giant, 
despite it developing a strong mean molecular weight 
gradient and forming a core-like structure (in terms 
of composition - the density remains somewhat lower 
than in the standard case). We therefore conclude that 
a mean molecular weight gradient alone is not suffi- 
cient to determine whether a star becomes a red giant 
or not, but it is necessary. 

So what other conditions must be met if we are to 
get a red giant? One obvious candidate is the rate of 
energy production. Our energy decoupled model shows 
that we can get redward motion if sufficient energy is 
deposited into the star. This should seem reasonable: 
in order to expand a gas, work has to be done and 
hence an energy source is needed. 

To investigate the role that energy generation may 
play in whether a star becomes a red giant or not, 
we make the following tests. Firstly, we start from a 
model which has a well-defined mean molecular weight 
gradient. We have chosen to use a model from the 
forced profile sequence, taken at a point when the mass 
of the hydrogen-exhausted core is 0.2 M©. We also 
stop hydrogen from being converted into helium, as 
we are trying to determine how the model reacts to 
changes in the energy generation only. Any changes in 
the chemical structure could potentially be confusing. 
Because of this, the models we obtain should not be 
regarded as 'evolutionary' models - they are merely 
the route the star takes in relaxing from one set of 
conditions to another. With these modifications, we 
are now free to alter the energy generation profile from 
that used in the forced profile sequence. 

The first change we make is to simply increase the 
energy generation profile by a fixed amount by mul- 
tiplying it by some factor. We also do the same to 
the neutrino loss rate profile. If we multiply these pro- 
files of a factor of between 2 and 5, the star contracts 
and moves blueward. If the energy generation profile 
is multiplied by a factor of ten, we obtain different be- 
haviour. The star expands and moves redward. How- 
ever, all the layers of the star are expanding so we do 
not get the desired compact core/extended envelope 
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Figure 8: Top panel: HR diagram showing the 
track of the model where energy is injected at 
0.3 M© (solid line), with the evolution being red- 
ward. The standard track is shown for comparison 
(dashed line). The crosses denote points at which 
the detailed profiles in the lower panel have been 
taken. Bottom panel: Radius as a function of mass 
for selected models in the sequence. Note that sig- 
nificant expansion occurs only above 0.2 M@. 

structure of a giant. 

One may then ask whether the way in which we 
deposit energy into the star is important. In a real red 
giant, energy is released in a shell rather than through- 
out the star as the above models have assumed. We 
therefore attempt the following. We take the energy 
generation profile from the 'forced profile' sequence 
and the same input model as above. This time, instead 
of multiplying the whole energy generation profile, we 
increase the energy generation rate in a well defined re- 
gion. This is done by adding a narrow Gaussian term 
at a given mass. We have tried placing the term at 
0.1, 0.3 and 0.5 Mq, which correspond to additional 
energy generation inside the hydrogen-exhausted re- 
gion (i.e. the 'shell') of the initial model, just outside 
the H-exhausted region, and within the envelope rep- 
sectively. 

The results for the case of energy injection at 0.3 Mq 
are shown in Figure [H] We plot the evolution of the 
radius profile for several consecutive models. We see 
that above 0.2 Mq there is a substantial increase in the 
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radius for a give mass. Crucially, we have expansion 
inside the region at which energy is injected and it is 
only within the H-exhausted region that little change 
in radius takes place. We obtain similar results if the 
energy is injected at 0.5 Mq. 

As discussed above, it appears that a mean molecu- 
lar weight difference is necessary for a star to become a 
giant and we can illustrate this with the following test. 
We select a model without a strong molecular weight 
gradient from the forced model sequence and inject 
energy using the same Gaussian profile as above. We 
have chosen a model about halfway through the main 
sequence. In this case, we find that the entire star ex- 
pands regardless of where we inject the energy. Thus 
we deem the change in mean molecular weight to be a 
necessary condition for forming a giant. 

3.1 A spanner in the works: mass 
dependence 

The model of llbenl (|l993h which has constant opac- 
ity and constant mean molecular weight becomes a gi- 
ant. This is clearly at odds with the picture presented 
above, in which the mean molecular weight gradient 
plays a key role. We note that Iben's model is for a 
star more massive than our model (his model is 5 M0) 
so it seems sensible to repeat our pseudohelium exper- 
iment for this mass of star. It is here that we run into 
a problem. The evolution of this model is shown in 
Figure [3] along with a standarcQ 5 Mq model. Our 
5 Mq model, like that of Iben, does indeed become a 
giant even in the case where we keep its mean molec- 
ular weight constant (i.e. by setting the atomic mass 
of the pseudohelium to 1.5). There is only one conclu- 
sion we can drawn for this test - the mean molecular 
weight is not crucial to a star becoming a giant in all 
masses of star. However, we note that if we evolve a 
5Mq model in a homogenous manner (i.e. ensuring 
that the whole star is mixed as we did in section [272]) . 
it contracts just as the 1 Mq model did. 

One of our critical tests in the case of the 1 Mq 
model was the one in which we forced the energy gener- 
ation profile to remain as it was on the main sequence 
(the 'forced profile' sequence of section I2.4.2p . This 
demonstrated that the mean molecular weight gradi- 
ent alone was not enough to produce a giant. We have 
re-run this test for the 5 Mq model, taking the zero- 
age main sequence energy generation profile for this 
mass of star. We find that this star does not become 
a giant, even though we are able to grow an extremely 
large core (over 80% of the star's mass) . It seems that 
an increase in the stars luminosity is indeed a neces- 
sary condition for a star to become a red giant. This 
is perhaps an unsuprising result: in order to expand 
something, work must be done and this requires an 
input of energy. 

We are still in the process of investigating the 5 Mq 
model in order to determine why it becomes a giant. 

3 That is, one that also uses only CNO burning, has no 
convective mixing and which burns hydrogen to normal he- 
lium in the same way that our standard 1 Mq model did. 
This is not a normal 5 Mq model! 
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Figure 9: HR diagram showing the evolution of a 
5 Mq model evolved with a pseudohelium atomic 
mass of 1.5 (dashed line). A standard 5 Mq model 
is shown for comparison (solid). Both models use 
only CNO burning reactions and no mixing has 
been permitted in convective regions. The stan- 
dard evolution has been terminated on the giant 
branch. 



However, the fact that a more massive star can be- 
come a giant when it has a constant mean molecular 
weight is extremely important. Most work on why 
stars become red giants is done under one fundamen- 
tal assumption, namely that there is only one way to 
make a giantQ and it is common to all stars, regard- 
less of mass (and perhaps other properties too, such 
as metallicity) . As such, many studies have only ex- 
amined one part icular mass of star (e.g. iRenzini et alJ 
119921 ; llbenl Il99l ) and one should therefore take their 
conclusions in that light. 

That different ma sses of star may beha ve differ- 
ently is hinted at by lEggleton et alJ ljl998h . Using 
poly tropic models, these authors demonstrated that 
a jump in the mean molecular weight could affect the 
Schonberg-Chandrasekher limit. They found that if 
the ratio of the core mean molecular weight to that 
of the envelope is less than 3 then the core mass can 
be arbitrarily large, though there is a limit to the core 
radius. We may speculate that our 1 and 5Mq mod- 
els behave in different ways because they sit on differ- 
ent sides of this (or a similar) threshold. The 1 Mq 
model may not become a giant if it does not have a 
mean molecular weight gradient because it is unable 
to form a core. This connection will be looked at in 
future work. In addition, the tests we have carried 
out should be repeated on a range of masses to see 
if there are further difference we have not yet noticed 
and which might yield clues as to why stars become 
red giants. 



4 The assumption may be correct: it may appear that 
there is more than one way to make a giant simply because 
we are ignorant of the factors that cause it. 
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4 Conclusions 

We have experimented with the inputs to a stellar 
model in order to work out why stars become red gi- 
ants. Our models show two things seem to be impor- 
tant in our 1 M© star: 

• A star must have a mean molecular weight gra- 
dient if it is to become a giant, but this is not a 
sufficient condition for it to do so. 

• Sufficient energy must be supplied in order for 
the star to become a giant. 

The mean molecular weight is merely a tool - the star 
must have enough power to make use of that tool. 

However, this cannot be the whole story as to why 
stars become red giants. Running the same battery of 
tests on a 5 M Q model yields a different result. Our 

5 M Q star can become a giant without a mean molec- 
ular weight gradient. However, sufficient energy must 
still be supplied to the star in order for it to expand 
just as in the case of our 1 Mq model. We suggest that 
there may be more than one way of making a giant and 
we will address this problem in future work. 
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